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This  paper  investigates  the  relationship  between  industrialization,  electricity  consumption  and  C02 
emissions  in  case  of  Bangladesh  using  quarter  frequency  data  over  the  period  of  1975-2010.  The  ARDL 
bounds  testing  approach  is  applied  to  examine  cointegration  in  the  presence  of  structural  breaks 
stemming  in  the  series.  The  causal  relationship  among  the  variables  is  explored  by  applying  the 
innovative  accounting  approach  (1AA). 

Our  results  indicate  that  the  variables  are  cointegrated  for  a  long  run  relationship.  We  find  that 
financial  development  adds  in  energy  pollutants.  Electricity  consumption  contributes  to  C02  emissions. 
Trade  openness  also  has  a  positive  impact  on  energy  pollutants.  The  results  unveil  that  EKC  is  existed 
between  industrial  development  and  C02  emissions  in  case  of  Bangladesh.  Our  causality  analysis  shows 
that  electricity  consumption  Granger  causes  energy  pollutants,  industrial  growth  and  financial  develop¬ 
ment.  The  unidirectional  causality  exists  running  from  financial  development  to  trade  openness  and 
trade  openness  Granger  causes  industrial  development.  This  study  opens  up  new  insights  for  policy 
makers  in  formulating  a  comprehensive  economic,  financial  and  trade  policy  to  sustain  industrialization 
by  improving  the  environmental  quality. 
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1.  Introduction 

In  recent  years,  rapid  industrialization,  increased  population 
and  significant  change  in  pattern  of  trade  and  financial  sector 
development,  the  threat  of  global  warming  and  climate  change  is 
increasing  for  the  last  two  decades  in  Bangladesh.  The  climate 
change  issue  is  fundamentally  important  for  Bangladesh.  The  EKC 
hypothesis  explains  that  the  level  of  carbon  emissions  is  expected 
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to  rise  for  many  of  the  economies,  contributing  further  to  global 
warming.  The  existence  of  the  EKC  is  a  debatable  issue  and 
depending  on  the  level  of  development  in  the  area  of  industria¬ 
lization,  energy  consumption,  and  pattern  of  trade  with  major 
trading  partners  might  differ  significantly  in  terms  of  their  output- 
emissions  nexus.  These  fundamentals  factors  are  associated  with 
emissions,  which  are  also  closely  related  to  the  relations  between 
energy  utilization  in  Bangladesh.  In  addition,  several  other  factors, 
e.g.  the  composition  of  growth,  type  of  economic  activities,  the 
industrialization  intensity  of  foreign  trade  and  financial  develop¬ 
ment  influence  the  emissions-energy-output  nexus  in  Bangladesh. 

In  Bangladesh  a  South  Asian  country,  a  number  of  industrial 
plants  were  established  in  different  areas  (Dhaka,  Chittagong, 
Narangong,  Gazipur,  etc.)  of  the  country  following  the  both  trade 
and  financial  liberalization  in  the  1990s;  however,  environmental 
issues  were  not  considered  seriously  at  that  time.  Thus,  the  natural 
environment  has  been  continuously  degrading  due  to  industrial 
enterprises'  unsustainable  practices  (raw  leather  processing  and 
chemical  wastage  and  deforestation  due  to  urbanization).  Accord¬ 
ing  to  the  World  Bank  [1]  information  for  Bangladesh,  it  has 
emitted  about  one  tenth  of  the  world's  C02  emissions  in  2006, 
despite  the  fact  that  its  160  million  people  represent  about  2.4%  of 
the  world's  population.  The  recent  threat  in  climate  change  for 
economic  development  strategies  in  many  developing  or  newly 
industrialized  countries  in  particular  for  poor  industrialization. 
Previous  research  study  by,  Rahman's  [2]  revealed  that  many 
people  are  being  affected  by  the  resultant  industrial  pollution  in 
Bangladesh.  Another  study  on  the  similar  sector  explains  the 
national  pollution  profile  also  shows  that  industrial  sub-sectors 
contribute  to  a  large  percentage  of  the  pollution  load  in  the 
country  [3],  Considering  the  information  reported  by  Bala  and 
Yusuf  [4],  suggested  that  Bangladesh's  existing  industrial  concerns, 
which  are  responsible  for  environmental  pollution  and  public 
health  hazards,  must  take  appropriate  measures  for  preventing 
its  pollution. 

In  this  study,  we  contribute  the  causal  relationship  between 
industrialization,  electricity  consumption  and  energy  pollutants  in 
case  of  Bangladesh  using  quarter  frequency  data  over  the  period  of 
1975-2010.  The  intensity  in  carbon  emissions,  the  demand  for 
electricity  consumption  and  rapid  industrialization  have  raised  some 
questions  in  Bangladesh's  perspective  such  as  (1)  Is  there  any  long¬ 
term  and  short-run  relationship  between  these  variables  in  Bangla¬ 
desh?  (2)  What  are  the  policy  implications  of  the  findings  for 
Bangladesh  and  other's  developing  or  newly  industrialized  econom¬ 
ics?  In  addition,  we  also  consider  the  influence  of  trade  and  finance 
within  the  emissions-energy-industrialization  nexus.  The  methodolo¬ 
gical  discussion  is  based  on  the  application  of  the  ARDL  bounds  testing 
approach  in  the  presence  of  structural  breaks.  These  methods  are 
performing  better  to  other,  especially  when  the  sample  size  is  small.  In 
addition,  the  causal  relationship  between  the  variables  is  explored  by 
applying  the  innovative  accounting  approach  (IAA).  Due  to  our 
imperfect  knowledge,  this  study  may  be  a  comprehensive  effort  on 
this  topic  for  the  economy  of  Bangladesh  and  it  has  five  ways 
contribution  to  the  energy  literature  by  applying:  (i)  Zivot-Andrews 
[5]  structural  break  unit  root  test;  (ii)  the  ARDL  bounds  testing 
approach  to  cointegration  for  long  run  relationship  between  the 
variables  in  the  presence  of  structural  breaks;  (iii)  OLS  and  ECM  for 
long  run  and  short  run  impacts;  (iv)  the  VECM  Granger  causality 
approach  for  causal  relationship;  and  (v)  innovative  accounting 
approach  (IAA)  to  test  the  robustness  of  causality  analysis.  The  main 
empirical  findings  of  this  study  found  that  a  long  run  relation  exists 
among  the  variables.  The  interesting  finding  is  that  the  EKC  is  existed 
between  industrial  development  and  C02  emissions  in  Bangladesh. 

The  remainder  of  this  paper  is  structured  as  follows:  section  2 
explains  literature  review.  In  Section  3,  we  outline  the  econometric 
specification  and  estimation  methodology  and  discuss  how  various 


hypotheses  are  tested,  while  Section  4  provides  a  discussion  of  our 
empirical  results.  Finally,  Section  5  discusses  major  findings  and 
concludes  the  paper. 

2.  Literature  review 

Research  strands  in  the  earlier  empirical  literature  to  examine 
the  relations  between  economic  growth,  energy  utilization  and 
emissions  [6],  The  previous  research  on  energy  economics  mainly 
concentrated  on  testing  of  the  existence  of  environmental  Kuznets 
curve  (EKC).  In  this  context,  the  pioneering  study  of  Kuznets  [7] 
which  claimed  for  an  inverted  U-shaped  relationship  between 
economic  growth  and  income  inequality  has  been  later  reformu¬ 
lated  to  test  similar  inverted-U  relationship  between  economic 
growth/income  and  environmental  quality.2 3  The  results  of  such 
research  are  however  contradictory  and  in  many  cases  researchers 
failed  to  establish  the  inverted-U  relationship.2 

The  second  argument  of  the  empirical  evidence  looks  at  the 
attention  between  energy  consumption  and  output  suggesting 
that  energy  consumption  and  income  may  be  jointly  determined. 
Following  the  seminal  work  of  Kraft  and  Kraft  [29];  several  others 
including  Masih  and  Masih  [30],  Yang  [31],  Wolde-Rufael  [32], 
Narayan  and  Singh  [33],  Narayan  et  al.[34]  tested  the  energy 
consumption  and  economic  growth  nexus  with  a  variety  of 
econometrics  techniques  and  for  different  panel  of  countries.  The 
existing  studies  on  energy-output  nexus  fall  into  four  broad 
categories:  For  instance,  studies  that  find  no  causality  between 
energy  consumption  and  income  in  case  of  a  U.S.  include  those  by 
Akarca  and  Long  [35],  Yu  and  Hwang  [36],  Yu  and  Jin  [37],  Stern 
[38],  and  Cheng  [39],  Fatai  et  al.  [40]  find  no  causality  between 
these  variables  in  New  Zealand.  The  same  result  is  true  in  the  UK 
[41]  and  France  [42],  In  contrast,  the  results  of  both  bidirectional 
and  unidirectional  causality  are  prevalent  particularly  in  develop¬ 
ing  nations  and  numerous  studies  have  found  such  bidirectional 
causality  with  real  life  data.4  Studies  find  two-way  causation  in 
South  Korea  include  those  of  Masih  and  Masih  [30],  Glasure  [43] 
and  Oh  and  Lee  [44],  A  similar  causation  for  the  Philippines  and 
Thailand  is  found  by  Asafu-Adjaye  [45].  Unidirectional  causation 
from  output  to  energy  consumption  is  found  in  different  coun¬ 
tries.5 6  A  number  of  studies  have  found  the  direction  of  causality 
from  energy  consumption  to  output  growth. 

Previous  researchers  have  attention  not  only  at  income  or 
economic  development  variables  but  also  extended  their  analysis 
to  include  other  aspects  such  as  financial  development  or  trade 
openness  or  trade  intensity  or  industry  value-added  of  a  country. 
These  factors  are  important  for  industrialized  countries  like 
Bangladesh.  The  relations  between  carbon  emissions  and  trade 


2  For  example,  Grossman  and  Krueger  [8],  Shafik  [9],  Dinda  and  Coondoo  [10], 
Heil  and  Selden  [11],  and  Friedl  and  Getzner  [12]  attempted  to  test  the  existence  of 
EKC  for  different  countries. 

3  The  studies  such  as  Ang  [13]  for  France;  Ang  [14]  and  Saboori  et  al.[15]  in 
Malaysia;  Dhakal  [16],  Tiwari  [17],  and  Tiwari  et  al.[  18]  in  India;  Shiyi  [19],  Jalil  and 
Mahmud  [20],  and  Zhang  and  Cheng  [21]  for  China;  Akbostanci  et  al.[22]  and 
Ozturk  et  al.[23]  in  Turkey;  Lean  and  Smyth  [24]  for  ASEAN  countries;  Narayan  and 
Narayan  [25]  for  43  developing  countries,  Fodha  and  Zaghdoud  [26]  for  Tunisia, 
Nasir  and  Rehrnan  [27]  and  Shahbaz  et  al.  [28(a)]  in  Pakistan;  Shahbaz  et  al.  [63]  for 
Romania  and,  Shahbaz  and  Leitao  [91  ]  for  Portugal  noted  the  empirical  evidence  of 
the  EKC. 

4  Some  of  the  studies  in  this  context  are  the  work  of  Shahbaz  and  Lean  [92]  and 
Shahbaz  et  al.[28(b)J  for  Pakistan;  Soytas  and  Sari  [49]  for  Argentina;  Ghali  and  El- 
Sakka  [93]  for  Canada;  Zachariadis  and  Pashourtidou  [94]  for  Cyprus;  Hondroyian- 
nis  et  al.  [95]  for  Greece;  Jumbe  [96]  for  Malawi;  Masih  and  Masih  [52]  for 
Pakistan,  etc. 

5  Like  Australia  [46],  India  [47],  Singapore  [48,49],  and  Taiwan  [50], 

6  These  studies  include  those  for  China  [51],  Fiji  [33],  India  [52],  Indonesia  [45], 
Philippines  [41],  Sri  Lanka  [53],  Taiwan  [31],  Turkey  [49, [54],  Venezuela  [55], 
Tunisia  [56],  and  Pakistan  [57]. 
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consider  fact  that  developed  economies  would  specialize  in 
human  or  physical  capital  intensive  activities  which  are  less 
emissions  more  intensive  than  developing  countries.  Trade  there¬ 
fore  may  result  in  increased  pollution  in  developing  countries  due 
to  increased  industrialization  of  these  emission  intensive  goods  in 
these  countries.  The  study  of  Grossman  and  Krueger  [8]  is 
pioneering  effort  in  this  regard  while  similar  research  question 
has  also  been  addressed  by  Lucas  et  al.  [58],  Wyckoff  and  Roop 
[59],  Nohrnan  and  Antrobus  [60],  etc.  The  earlier  findings  fail  to 
conclude  the  conclusive  evidence  of  the  relationships  between 
trade  and  environmental  quality.  Applying  Turkish  data,  Halicioglu 
[61]  estimates  the  causal  relationship  between  trade,  C02  emis¬ 
sion,  and  income  and  energy  consumption  and  found  that  income 
is  the  most  crucial  determinant  of  carbon  emissions,  followed  by 
energy  consumption  and  finally  trade.  Using  the  similar  approach, 
the  importance  of  foreign  trade  in  determining  the  level  of  C02 
emissions  has  also  been  emphasized  by  Anderson  et  al.  [62],  In 
their  analysis,  the  authors  find  that  the  exports  of  China  played  an 
important  role  in  generating  emissions  in  the  transport  sector 
which  is  greater  than  the  emissions  attributable  to  imports.  Recent 
study  by  Shahbaz  et  al.  [28(b)]  also  found  that  trade  openness 
increases  C02  emissions  due  to  poor  implementation  of  environ¬ 
mental  regulation  in  case  of  Pakistan.  Furthermore,  financial 
instability  is  also  positively  linked  with  environmental  degrada¬ 
tion  [63], 

Applying  augmented  vector  autoregression  (VAR)  approach  of 
the  Toda  and  Yamamoto  [64]  after  incorporating  gross  fixed  capital 
formation  and  labor  force  into  the  model,  Soytas  et  al.  [6] 
investigated  energy  consumption,  output  and  carbon  emissions 
nexus  for  the  USA  and  found  no  causal  relationship  between 
income  and  carbon  emissions  or  between  energy  use  and  income. 
Using  the  similar  procedure,  Zhang  and  Cheng  [21]  investigated 
the  energy  consumption,  output  and  carbon  emissions  nexus  for 
China,  controlling  for  capital  and  urban  population.  The  authors 
found  unidirectional  long-run  causality  running  from  GDP  to 
energy  consumption  and  from  energy  consumption  to  carbon 
emissions.  Considering  the  study  by  Sari  and  Soytas  [65],  investi¬ 
gated  the  relationship  between  carbon  emissions,  income,  energy 
and  total  employment  in  five  OPEC  countries,  the  authors  found 
that  none  of  the  following  countries  namely  Algeria,  Indonesia, 
Nigeria,  Saudi  Arabia  and  Venezuela  need  to  sacrifice  economic 
growth  in  order  to  reduce  C02  emissions.  Applying  the  ARDL 
cointegration  approach;  Halicioglu  [61]  examined  the  cointegra¬ 
tion  in  a  log  linear  quadratic  relationship  between  per  capita  C02 
emissions,  per  capita  energy  use,  per  capita  real  income,  square  of 
per  capita  real  income  and  openness  ratio,  and  reported  that  there 
is  a  short-and-long  runs  bidirectional  causality  between  carbon 
emissions  and  income  in  Turkey.  In  a  similar  study;  Jalil  and 
Mahmud  [20]  found  unidirectional  causality  from  economic 
growth  to  C02  emissions  in  China.  The  study  also  indicates  that 
carbon  emissions  are  mainly  determined  by  income  and  energy 
consumption  in  the  long-run.  Trade  has  a  positive  but  statistically 
insignificant  impact  on  C02  emissions. 

Using  variables  like  stock  market  value  added,  foreign 
direct  investment,  deposit  money  bank  assets,  capital  account 
convertibility,  financial  liberalization,  financial  openness  etc,  in 
order  to  capture  the  level  of  financial  development,  Tamazian  et  al. 
[66]  conducted  detailed  analysis  on  the  relationship  between 
financial  development  and  environmental  degradation.  Over  time 
the  research  on  environment-development  nexus  has  not  only 
extended  and  modified  in  terms  of  the  research  questions 
addressed,  but  also  has  developed  from  the  viewpoint  of  meth¬ 
odology  and  econometric  modeling  techniques.  Most  of  the 
studies  like  that  of  Peng  and  Sun  [67]  and  Akbostanci  et  al.  [22] 
have  applied  cointegration  technique  and  Granger  causality 
method  to  understand  the  nexus.  In  terms  of  the  findings  of 


these  analyses,  no  clear  cut  conclusion  can  be  made  and  results 
differ  depending  upon  the  variables  used  and  countries  consid¬ 
ered.  In  case  of  ASEAN  countries,  a  recent  study  conducted 
by  Lean  and  Smyth  [24]  with  panel  data  attempted  to  examine 
the  long-run  relationship  between  C02  emission,  electricity 
consumption  and  output  as  well  as  the  causal  relationship 
between  these  variables.  Their  analysis  revealed  that  there  is 
short-run  panel  Granger  causality  from  carbon  emissions  to 
electricity  consumption  and  long-run  unidirectional  causality 
from  electricity  consumption  and  C02  emissions  to  GDP.  Akbos¬ 
tanci  et  al.  [22]  also  found  pollution  and  income  variables  being 
cointegrated. 

Previous  study  on  Bangladesh's  energy  and  output  by  Mozum- 
der  and  Marathe  [68],  over  the  1971-1999  period,  found  that  per 
capita  GDP  causes  per  capita  electricity  consumption,  but  the 
reverse  is  not  true.  However,  small  sample  size  is  the  limitation  of 
this  study.  From  the  annual  data  from  1971  to  2007,  Uddin  et  al. 
[69]  investigated  the  inter-temporal  causal  relationship  between 
energy  consumption  and  economic  growth  in  Bangladesh  and 
unidirectional  causality  runs  from  energy  consumption  to  eco¬ 
nomic  growth  in  Bangladesh  and  then  restrictions  on  the  use  of 
energy  could  lead  to  a  reduction  in  economic  growth.  The  limita¬ 
tion  of  this  approach  is  bivariate  framework  to  investigate  the 
relationship  between  energy  and  income.  A  study  in  Bangladesh 
by  Paul  and  Uddin  [70]  examined  energy-output  dynamics  for 
Bangladesh  and  found  that  while  fluctuations  in  energy  consump¬ 
tion  do  not  affect  output  fluctuations,  movements  in  output 
inversely  affect  movements  in  energy  use.  The  results  of  Granger 
causality  tests  in  this  respect  are  consistent  with  those  of  innova¬ 
tive  accounting  that  includes  variance  decompositions  approach 
and  impulse  responses.  Autoregressive  distributed  lag  models  also 
suggest  a  role  of  output  in  Bangladesh's  energy  use.  However,  Paul 
and  Uddin  study  does  not  consider  the  effect  of  openness,  finance 
and  industrialization  impact  in  their  assigned  model  specification. 
Our  study  filled  this  gap.  Working  with  annual  data,  Ahamad  and 
Islam  [71]  examine  the  causal  relationship  between  per  capita 
electricity  consumption  and  per  capita  GDP  of  Bangladesh  using 
the  vector  error  correction  specified  Granger  causality  for  the 
period  of  1971-2008.  Empirical  findings  reveal  that  there  is  short- 
run  unidirectional  causal  flow  running  from  per  capita  electricity 
consumption  to  per  capita  GDP  without  feedback.  Structural  break 
was  absent  in  their  estimation  process  which  is  valid  for  Bangla¬ 
desh  energy  output  dynamics  [70],  Ahiduzzaman  and  Islam  [72] 
study  presents  a  review  of  the  potential  and  utilization  of  the 
renewable  energy  sources  in  Bangladesh.  In  case  of  Bangladesh, 
Alam  et  al.  [73]  investigated  the  relationship  between  energy 
consumption,  C02  emissions  and  economic  growth  by  applying 
Johansen  and  Juselius  [74],  cointegration  and  the  VECM  Granger 
causality  approaches.  Their  results  noted  that  energy  consump¬ 
tion  Granger  causes  economic  growth  and  C02  emissions. 
Furthermore,  economic  growth  is  Granger  cause  of  C02  emis¬ 
sions.  Amin  et  al.  [75],  however,  found  a  neutral  effect  between 
C02  emissions  and  economic  growth.  Both  of  these  earlier 
papers  have  shortcomings  in  the  methodological  framework 
that  structural  break  is  not  considered  in  their  dynamics 
estimation  which  is  filled  up  in  this  current  study.  Recent 
studies,  the  top-ten  environment  polluting  industries  in  Bangla¬ 
desh;  these  are  the  tannery,  pulp  and  paper,  fertilizer,  textile 
and  cement  industries,  Hoque  and  Clark  [76]  using  both  primary 
and  secondary  data  from  each  industry  group,  two  sample 
plants  were  selected  with  five  executives  participating  from 
each  plant.  This  study  highlights  the  reality  of  Bangladeshi 
industrial  plants  in  applying  pollution  prevention  initiatives.  It 
reveals  that  compared  to  leading  firms  in  developed  countries, 
pollution  prevention  initiatives  in  Bangladesh  are  considerably 
underutilized. 
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3.  Data  construction  and  methodology 

3.1.  Data  construction 

The  data  used  for  the  analysis  has  mainly  been  taken  from 
World  Development  Indicator  (CD-ROM,  2012),  published  by  the 
World  Bank.  It  covers  the  period  of  1975-2010.  The  variables 
under  consideration  are  carbon  emissions,  industrial  value-added 
as  a  share  of  GDP,  electricity  consumption,  domestic  credit  to 
private  sector  as  a  share  of  GDP  and  trade  openness  (exports + 
imports)  as  a  share  of  GDP.  For  our  analysis  these  variables  have 
been  converted  to  logarithmic  form.  Here,  carbon  emissions  is 
measured  in  per  capita,  electricity  use  is  calculated  as  Electric 
power  consumption  (kWh)  per  capita,  industrial  value  added 
measure  in  real  per  capita,  financial  development  is  measured  in 
domestic  credit  to  the  private  sector  per  capita  in  real  terms.  As  a 
proxy  for  foreign  trade,  trade  openness  is  calculated  as  the  sum  of 
real  export  and  imports  per  capita  has  been  used.  Thus,  we  get  140 
observations  on  each  series  ranging  from  1975  to  2010,  the  longest 
possible  joint  dataset  on  Bangladesh.  We  have  converted  all  the 
annual  series  into  quarterly  data  to  avoid  the  problem  of  degree  of 
freedom  and  efficient  empirical  results.  We  used  quadratic  match 
sum  method  to  transform  all  the  variables  into  quarter  frequency 
following  Romero  [77]  and  McDermott  and  McMenamin  [78].'7 

3.2.  Zivot-Andrews  unit  root  test  logarithm 

Zivot-Andrews  [5]  developed  a  unit  root  test  in  the  presence  of 
structural  break  in  the  econometric  models.  The  model  is  applic¬ 
able  in  the  context  of  macroeconomic  series.  The  model  con¬ 
structed  at  the  three  different  forms  such  as  level,  function  and 
both  in  trend  and  intercept  function  of  the  variables  to  be  applied 
in  the  time  series  analysis. 

3.3.  The  ARDL  bounds  testing  approach  to  cointegration 

The  empirical  analysis  of  this  study  employs  the  ARDL  bounds 
testing  approach  [79,  80]  to  explore  the  existence  of  long-run 
equilibrium  relations  among  the  variables.  This  approach  has 
number  of  advantages  compared  to  other  traditional  cointegration 
techniques.  First,  it  allows  for  smaller  sample  sizes.  Second,  it  can 
be  used  regardless  of  whether  the  variables  are  purely  I  (0),  purely 
I  (1),  or  mutually  cointegrated.  Third,  it  provides  unbiased  long- 
run  estimates  and  valid  t-statistics.  Finally,  this  approach  provides 
a  method  of  assessing  short-run  and  long-run.  Moreover,  a 
dynamic  unrestricted  error  correction  model  (UECM)  can  be 
derived  through  a  simple  linear  transformation.  The  UECM  inte¬ 
grates  the  short-run  dynamics  with  the  long-run  equilibrium 
without  losing  any  long-run  information.  The  UECM  is  specified 
as  follows: 

A  In  Ct  =  a !  +aDUMDUM+ac  In  Ct-\+(Xi  In It^i+aE  In Ec_ r 

p 

+aF  lnFt_i+a0  In  Ot_i+  £  atA  In  Ct_; 

i  =  i 

+  £  cijA  In  It~j+  2  akAh\E,_k+  £  atA  In  Ft_t 
]=o  k  = o  l=o 

t 

+  2  amA  In  Or_m+/rf  (1) 

m  =  0 

where  In  Ct  =  natural  log  of  C02  emissions  per  capita;  In  It  =  natural 
log  of  real  industrial  value-added  per  capita;ln  Et  =  natural  log  of 


7  All  the  series  are  transformed  into  natural-log  form.  The  log-linear  specifica¬ 
tion  of  the  series  provides  efficient  and  consistent  empirical  results  compared  to 
simple  specification.  Furthermore,  transformation  of  series  into  natural  log  reduces 
the  sharpness  of  the  data  and  provides  reliable  empirical  findings. 


electricity  consumption  per  capita;  In  Ft  =  natural  log  of  real  domestic 
credit  to  private  sector  per  capita;  In  Ot  =  natural  log  of  real  trade 
openness  ( exports + imports)  per  capita  and  DUM  is  dummy  variable 
to  capture  structural  break  based  on  Z-A  unit  root  test.  A  is  the  first 
difference  operator,  T  is  the  time  trend  and  /qis  the  error  term.  The 
optimal  lag  structure  of  the  first  difference  regression  is  selected  based 
on  Akaike  Information  Criteria  (AIC).  The  F-test  is  used  in  bounds  test 
for  the  existence  of  the  long-run  relationship  [80]  and  it  tests  for  the 
joint  significance  of  lagged  level  variables  involved.  The  null  hypoth¬ 
esis  of  the  non-existence  of  a  long-run  relationship  for  the  equation  of 
being  H0  :  ac  =  cti  =  aE  =  aE  =  a0  =  0  against  the  alternative  hypoth¬ 
esis  Ha  :  ac  ¥=  ^  aE  =£  aF  a0  ^  0.  According  to  the  ARDL 

approach,  if  the  F-statistic  exceeds  the  upper  critical  value,  we 
conclude  in  the  favor  of  long  run  relationship.  If  the  F-statistic  falls 
below  the  lower  critical  value,  we  cannot  reject  the  null  hypothesis  of 
no  cointegration.  However,  if  the  F-statistic  lies  between  the  two 
bounds,  inference  is  inconclusive.8 

3.4.  Innovative  accounting  approach  for  Granger  causality 

According  to  the  Shahbaz  et  al.  [28]  the  Granger  causality  test  is 
unable  to  indicate  how  much  feedback  exists  from  one  variable  to 
the  other.  Due  to  the  limitation  of  the  Granger  causality  test,  we 
include  innovative  accounting  approach  (IAA)  to  investigate  the 
dynamic  causal  relationships  among  carbon  emissions,  electricity 
consumption  and  industrialization,  financial  development  and 
trade  openness.  The  uniqueness  of  the  IAA  is  that  it  avoids  the 
problem  of  endogeneity  and  integration  of  the  series.  This 
approach  has  an  advantage  compared  to  the  VECM  Granger  cau¬ 
sality  test  because  the  latter  only  shows  a  causal  relationship 
between  the  variables  within  the  sample  period  while  the  former 
illustrates  the  extent  of  causal  relationship  ahead  the  selected 
sample  period.  The  IAA  includes  variance  decomposition  method 
and  impulse  response  function.  This  procedure  decomposes  fore¬ 
cast  error  variance  for  each  series  following  a  standard  deviation 
shock  to  a  specific  variable  and  enables  us  to  test  which  series  is 
strongly  impacted  and  vice  versa. 

For  instance,  if  a  shock  in  C02  emissions  has  significant  effects 
of  industrialization  but  a  shock  occurring  in  industrialization  only 
affect  very  minimum  the  variations  of  C02  emissions.  Then,  this  is 
inferred  as  a  unidirectional  causality  runs  from  C02  emissions  to 
industrialization.  If  the  industrialization  explains  more  of  the 
forecast  error  variance  of  C02  emissions;  then  we  deduce  that 
industrialization  causes  C02  emissions.  The  bidirectional  causality 
exists  when  shocks  in  C02  emissions  and  industrialization  have  a 
strong  impact  on  the  variability  of  C02  emissions  and  industriali¬ 
zation  respectively.  If  shocks  occur  in  both  series  do  not  have  any 
impact  on  the  industrialization  and  C02  emissions  then  there  is  no 
causality  between  the  variables. 

Impulse  response  function  helps  us  to  trace  out  the  time  path 
of  the  impacts  of  shocks  of  variables  in  the  VAR.  One  can 
determine  how  much  C02  emissions  responses  due  to  its  own 
shock  and  shock  in  industrialization.  We  support  the  hypothesis 
that  industrialization  causes  C02  emissions  of  the  impulse 
response  function  indicates  significant  response  of  C02  emissions 
to  shocks  in  industrialization  than  other  variables.  A  strong  and 
significant  reaction  of  C02  emissions  to  shocks  in  industrialization 
implies  that  C02  emissions  cause  industrialization.  This  study 
incorporates  C02  emissions,  industrialization,  electricity  consump¬ 
tion,  financial  development  and  trade  openness  to  examine  the 
relationship  between  C02  emissions  and  its  determinants  in  the 


8  In  addition,  Pesaran  et  al.[80]  caution  that  critical  values  for  the  bound  test 
are  sensitive  to  the  number  of  regressors  (k)  in  the  model,  and  Narayan[81 1  argues 
that  critical  values  of  F-test  depend  on  sample  size. 
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VAR  model.  AVAR  system  takes  the  following  form  [82]: 

Vt=  £  5iVt_,+ifc  (2) 

i  =  1 

where 

Vt  =  (Ct,It,Et,Ft,  Ot) 

Sj  are  the  estimated  coefficients  and  ;/  is  a  vector  of  error  terms. 

4.  Results  and  discussions 

Descriptive  statistics  and  correlation  matrix  is  presented  in 
Table  1.  Based  on  Jarque-Beratest  statistics,  our  results  indicate 
that  all  the  series  are  normally  distributed  having  zero  mean  while 
variance  is  constant.  This  implies  that  we  should  peruse  for  further 
analysis.  The  correlation  matrix  reveals  a  positive  association 


Table  1 

Descriptive  statistics  and  correlation  matrix. 


Variables 

In  Ct 

In  It 

In  Et 

In  Ft 

In  Of 

Mean 

-1.9418 

7.8960 

4.0189 

4.3215 

8.0852 

Std.  Dev. 

0.4897 

0.4643 

0.8585 

0.5115 

0.6062 

Jarque-Bera 

2.4773 

2.2765 

2.0711 

1.1488 

2.3841 

Probability 

0.2897 

0.3203 

0.3550 

0.5630 

0.3035 

In  C, 

1.0000 

In  It 

0.9832 

1.0000 

In  Et 

0.9898 

0.9915 

1.0000 

In  Ft 

0.9845 

0.9746 

0.9793 

1.0000 

In  Of 

0.9261 

0.9599 

0.9341 

0.9221 

1.0000 

between  the  underlying  variables.  For  instance,  industrialization 
is  positively  linked  to  C02  emissions  and  same  finding  is  for 
electricity  consumption  and  C02  emissions.  A  positive  correlation 
is  found  of  financial  development  and  trade  openness  with  C02 
emissions.  Electricity  consumption,  financial  development  and 
trade  openness  are  positively  related  with  industrialization.  A 
positive  correlation  of  financial  development  and  trade  openness 
is  found  with  electricity  consumption.  There  is  a  positive  correla¬ 
tion  between  financial  development  and  trade  openness. 

There  are  numerous  unit  root  tests  available  to  test  stationarity 
properties  of  time  series  data.  These  unit  root  tests  are  ADF,  PP,  DF- 
GLS,  KPSS  and  Ng-Perron.  Generally,  these  unit  root  tests  are  not 
reliable  in  the  presence  of  structural  break  [83],  This  limitation  of 
classical  unit  root  tests  such  as  ADF,  PP,  DF-GLS  and  KPSS  has  been 
covered  by  applying  Zivot-Andrews  [5]  structural  break  unit  root 
test.  Zivot-Andrews  unit  root  test  contains  information  about  one 
unknown  structural  break  in  the  series.  The  results  of  Zivot  and 
Andrew  [5]  unit  root  test  are  presented  in  Table  2.  This  empirical 
evidence  indicates  that  the  series  is  non-stationary  at  level  but 
found  to  be  stationary  at  1st  difference.  This  implies  that  all  the 
series  are  integrated  at  1(1). 

According  to  the  ARDL  approach,  lag  order  of  the  variables  is 
important  for  the  model  specification.  Table  3  indicates  the  lag 
length  criterion.  In  this  paper,  we  followed  Akaike  information 
criterion  to  select  an  appropriate  lag  length.  It  is  pointed  by 
Liitkepohl  [84]  that  A1C  has  superior  power  properties  for  small 
sample  data  compared  to  any  lag  length  criterion.  Akaike  informa¬ 
tion  criterion  provides  efficient  and  consistent  results  as  compared 
to  final  prediction  error  (FPE),  Schwarz  information  criterion  (SBC) 
and  Hannan-Quinn  information  criterion  (HQ).  Based  on  the 
empirical  evidence  provided  by  AIC,  we  find  that  the  optimum 
lag  is  6  in.  such  quarter  frequency  data  over  the  period  of  1975- 
2010  in  case  of  Bangladesh. 


Table  2 

Z-A  unit  root  analysis 


Variable 

Level 

1st  Difference 

T-statistic 

Time  break 

Decision 

T-statistic 

Time  break 

Decision 

In  C, 

-4.644  (2) 

1994Q1 

Unit  root 

-10.513  (3)* 

1984Q1 

Stationary 

In  It 

-4.940  (3) 

1979Q3 

Unit  root 

-15.552  (3)* 

1978Q4 

Stationary 

In  /[ 

-4.449  (3) 

1979Q3 

Unit  root 

15.422  (4)* 

1978Q4 

Stationary 

In  E, 

-4.420  (2) 

2004Q2 

Unit  root 

-9.120  (3)* 

1980Q3 

Stationary 

In  Ft 

-4.481  (2) 

1989Q2 

Unit  root 

-7.134  (4)* 

1979Q4 

Stationary 

In  Ot 

-4.142  (3) 

1985Q2 

Unit  root 

9.413  (5)* 

1984Q2 

Stationary 

*  Significant  at  1%  level  of  significance.  Lag  order  is  shown  in  parenthesis. 

Table  3 

Lag  order  selection  criteria. 

VAR  lag  order  selection  criteria 

Lag 

Log  L 

LR 

FPE 

AIC 

SC 

HQ 

0 

1760.621 

NA 

2.51  e- 19 

-25.80325 

-25.67475 

-25.75103 

1 

3399.985 

3109.969 

1.44e  — 29 

-49.38213 

-48.48263 

-49.01659 

2 

3552.110 

275.1679 

2.62e-30 

-51.08985 

-49.41936 

-50.41101 

3 

3567.684 

26.79737 

3.57e-30 

-  50.78948 

-48.34799 

-49.79732 

4 

3578.486 

17.63231 

5.24e-30 

-50.41892 

-47.20643 

-49.11344 

5 

3782.955 

315.7243 

4.51  e- 31 

-52.89640 

-48.91292 

-51.27762 

6 

3937.214 

224.5824* 

8.20e  — 32* 

-  54.63550* 

-49.88102* 

—  52.70340* 

7 

3957.061 

27.14392 

1.09e  — 31 

-  54.39796 

-48.87248 

-52.15254 

8 

3970.924 

17.73600 

1.62e  — 31 

-  54.07241 

-47.77593 

-51.51368 

Note:  LR,  sequential  modified  LR  test  statistic  (each  test  at  5%  level);  FPE,  final  prediction  error;  AIC,  Akaike  information  criterion;  SC,  Schwarz  information  criterion;  HQ, 
Hannan-Quinn  information  criterion. 

*  Lag  order  selected  by  the  criterion. 
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Table  4 

ARDL  cointegration  analysais 


Variable 

In  Ct 

In  It 

In  I2 

In  Et 

In  Ft 

In  TR, 

Structural  break 

1994Q1 

1979Q3 

1979Q3 

2004Q2 

1989Q2 

1985Q2 

F-statistics 

5.202* 

3.372*** 

3.805*** 

4.740** 

2.460 

6.487* 

Critical  values'1 

1  %  level 

5  %  level 

10  %  level 

Lower  bounds 

3.60 

2.87 

2.53 

Upper  bounds 

4.90 

4.00 

3.59 

Diagnostic  test 

R2 

0.6720 

0.9847 

0.9845 

0.7650 

0.8356 

0.9116 

Adj-R2 

0.3440 

0.9757 

0.9754 

0.6269 

0.7389 

0.8516 

F-statistic 

3.0498** 

11.1130* 

12.5765* 

5.5369* 

8.6434* 

11.5461* 

*  Significant  at  1%. 

**  Significant  at  5%. 

***  Significant  at  10%.level. 

a  Critical  values  bounds  are  from  [81]  with  unrestricted  intercept  and  unrestricted  trend. 


In  this  section  we  present  the  empirical  findings  based  on  our 
methodology  discussed  in  the  last  section.  The  estimation  results 
for  cointegration  are  presented  in  Table  4.  The  uniqueness  of  this 
paper  is  that  the  structural  break  incorporated  into  the  ARDL 
equation  while  computing  F-statistics.  If  the  calculated  F-statistic 
exceeds  the  upper  bound,  then  null  hypothesis  of  no  cointegration 
among  variables  in  In  Ct  can  be  rejected.  If  the  calculated 
F-statistic  falls  below  the  lower  bound,  the  null  hypothesis  of  no 
long-run  relationship  cannot  be  rejected  (see  footnote  3).  The 
results  reported  in  Table  4  show  that  there  is  evidence  of 
cointegration  when  In  It  and  In  Et  (In  Ot)  are  taken  as  dependent 
variables  in  the  presence  of  structural  break  at  1%,  10%  and  5%  (1%) 
except  In  Ft  in  the  model  in  case  for  Bangladesh.  This  shows  that 
there  are  five  cointegrating  vectors  validating  the  existence  of  long 
run  relationship  between  the  variables  in  presence  of  structural 
break  stemming  in  the  series  of  C02  emissions,  industrialization, 
electricity  consumption,  financial  development  and  trade  open¬ 
ness.  The  structural  breaks  in  1994Q1,  1979Q3,  2004Q2,  1989Q2 
and  1985Q2  reveal  that  Bangladesh  economy  underwent  the 
structural  adjustment  policies  in  the  early  1980s.  This  includes 
institutional  as  well  as  policy  level  alterations.  The  financial  and 
external  sectors  underwent  substantial  reforms  in  the  mid  1980s 
and  early  1990s  that  are  likely  to  have  influenced  real  output.  The 
real  output  is  also  likely  to  have  structural  breaks  as  a  consequence 
of  these  reforms.  In  the  past  decade  or  so,  Bangladesh  economy's 
growth  hovered  around  5-6%  per  annum.  From  the  restoration  of 
parliamentary  democracy  in  the  1990  and  onwards,  it  has  been 
observed  that  the  periodic  growth  averages  are  much  higher  than 
the  previous  two  decades.  GDP  growth  has  been  much  more 
volatile  till  1990,  which  became  much  more  stable  and  increasing 
in  the  subsequent  years. 

The  marginal  impacts  of  industrialization,  electricity  consump¬ 
tion,  financial  development  and  trade  openness  on  C02  emissions 
are  reported  in  Table  5.  The  results  indicated  that  all  the  coeffi¬ 
cients  are  according  to  our  expectations.  We  find  that  there  is  an 
inverted  U-shaped  relationship  found  between  industrialization 
and  carbon  emissions.  The  linear  and  squared  terms  of  industria¬ 
lization  have  positive  and  negative  relationships  with  carbon 
emissions.  We  note  that  a  1%  rise  in  real  industrial  value  added 
per  capita  will  raise  C02  emissions  by  4.0646%  while  a  negative 
sign  of  the  squared  term  of  real  industrial  value  added  per  capita 
seems  to  document  the  declining  of  energy  emissions  and  real 
industrial  value  added  per  capita  at  the  higher  level  of  income  is 
-0.2434%  at  the  1%  level  of  significance.  These  evidences  support 
the  EKC  hypothesis  revealing  that  carbon  emissions  increase  in  the 
initial  stage  of  industrialization  and  start  to  decline  once  industrial 
development  is  completed. 


Table  5 

Long  run  regression  analysis 


Dependent  variable  In  Ct 

Variable 

Coefficient 

Std.  Error 

t-  Statistic 

Constant 

-5.3609* 

0.6162 

-8.6988 

In  It 

4.0646* 

0.5883 

6.9087 

In  I2 

0.2.434 

0.0320 

-7.5881 

In  Ft 

0.1527* 

0.0550 

2.7752 

In  Et 

0.3127* 

0.0612 

5.1076 

In  Ot 

0.0893* 

0.0307 

2.9089 

R2 

0.9895 

Adj-R2 

0.9891 

F-  statistic 

26.0365* 

Diagnostic  test 

F-statistic 

P-value 

x2normal 

0.5147 

0.4624 

Z2 SERIAL 

0.4347 

0.5000 

X2ARCH 

0.3111 

0.7787 

X2  WHITE 

0.9534 

0.3230 

X2REMSAY 

0.1680 

0.8835 

Note:  x2  NORMAL  is  for  normality  test,  x2  SERIAL  for  LM  serial  correlation  test,  x2 
ARCH  for  autoregressive  conditional  heteroskedasticity,  x2  WHITE  for  white  hetero- 
skedasticity  and  x2  REMSAY  for  Resay  Reset  test. 

*  The  significant  at  1%  level. 

The  impact  of  electricity  consumption  of  C02  emissions  is  positive 
and  significant  at  1%  significance  level.  We  infer  that  all  else  is  same,  a 
1%  rise  in  electricity  consumption  is  linked  with  a  0.3127%  increase  in 
C02  emissions.  To  develop  low-carbon  society  in  Bangladesh,  a 
package  of  plans  of  Bangladesh,  2010-2021  such  as  National  Energy 
policy,  National  Renewable  Energy  policy  and  Strategic  Transportation 
Plans  was  introduced.  These  policies  are  emphasized  on  fuel  switch, 
introduction  of  renewable  and  nuclear  energy  in  power  sector, 
improvise  the  energy  efficient  equipment,  increment  of  public  trans¬ 
port  and  improve  traffic  management  system.  Realization  and  imple¬ 
mentation  of  such  policy  package  for  low-carbon  society  development 
requires  necessary  information  on  changes  in  socioeconomic  para¬ 
meters,  energy  demand,  C02  emissions  and  emissions  reduction  in 
future.  This  report  will  provide  essential  information  on  socioeco¬ 
nomic  parameters  and  potential  low-carbon  measures  in  making 
policy  plan  for  Bangladesh. 

Financial  development  and  C02  emissions  are  positively  linked 
and  statistically  significant  at  the  1%  level  of  significance.  It  is  noted 
that  a  0.15275%  increase  in  C02  emissions  will  be  increased  by  1%  rise 
in  financial  development.  Reforms  in  the  financial  sector  in  Bangla¬ 
desh  started  in  the  early  1980s  and  gained  the  pace  in  the  1990s.  The 
main  focus  of  these  reforms  was  to  improve  the  process  of  financial 
intermediation  by  taking  up  a  series  of  measures  related  to  legal, 
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policy  and  institutional  restructuring.  The  first  phase  of  reforms  in 
1980s  includes  denationalization  of  public  banks  in  1984,  allowing 
new  private  banks  in  1986,  establishment  of  a  National  Commission 
on  Money,  Banking  and  Credit  to  identify  problems  in  the  banking 
sectors  and  prescribe  policies  as  remedial  measures.  In  the  later 
phase  of  reforms,  the  government  allowed  for  market-determined 
deposits  and  lending  rates.  Other  measures  include  the  introduction 
of  indirect  monetary  instruments  to  replace  direct  credit  control, 
improvement  of  the  capital  base  of  commercial  banks,  and  reforms  in 
the  legal  framework  of  debt  recovery  [85],  In  1997,  the  ADB  approved 
a  program  loan  of  S80  million  that  was  aimed  at  enhancing  market 
capacity,  and  developing  a  fair,  transparent,  and  efficient  capital 
market  [86],  This  implies  that  environmental  quality  was  and  is  not 
on  the  priority  by  the  financial  sector  in  Bangladesh. 

Trade  openness  has  positive  and  statistically  significant  impact 
on  C02  emissions.  Keeping  other  economic  agents  constant,  a  1% 
increase  in  trade  openness  increases  C02  emissions  by  0.0893%.  By 
the  end  of  the  1970s,  Bangladesh  partially  changed  its  anti-exports 
bias  policies  and  by  the  mid-1980s  it  undertook  policies  and 
programs  that  resulted  in  consistent  improvement  in  the  incentive 
to  export.  By  the  1990s  Bangladesh  became  more  exports  oriented 
and  significant  improvements  have  been  made  in  exports  policy 
and  administration.  During  the  first  half  of  the  decade,  liberal¬ 
ization  policy  taken  by  the  governments  for  swelling  imports  gave 
support  to  the  flow  of  inputs  for  exports-oriented  and  domestic 
industries,  and  led  to  a  surge  in  consumer  and  non-production 
related  imports  [87],  This  shows  that  Bangladesh  followed 
exports-oriented  policies  at  the  cost  of  environment. 

Log  run  model  fulfills  the  assumptions  of  the  classical  linear 
regression  model  (CLRM)  such  as  normality  of  the  error  term, 
serial  correlation,  autoregressive  conditional  heteroskedasticity  as 
well  as  white  heteroskedasticity.  Our  results  find  that  residual 
term  shows  normal  distribution  and  evidence  of  no  serial  correla¬ 
tion  is  found.  The  autoregressive  conditional  heteroskedasticity 
and  white  heteroskedasticity  are  not  found.  Finally,  the  short  run 
model  is  well  specified  as  confirmed  by  Ramsey  RESET  tests. 

The  short  run  results  are  illustrated  in  Table  6.  The  results  find 
that  linear  and  nonlinear  terms  of  industrial  value-added  per 


Table  6 

Short  run  regression  analysis 


Dependent  variable: 

In  Cf 

Variable 

Coefficient 

Std.  error 

t-Statistic 

Constant 

0.0029* 

0.0006 

4.2912 

In  It 

0.3083 

1.8149 

0.1698 

In  I2 

-0.0235 

0.1201 

-0.1962 

In  Ft 

0.1265*** 

0.0760 

1.6629 

In  £t 

-0.1595** 

0.0673 

-2.3696 

In  Of 

0.1023** 

0.0403 

2.5342 

ECMt_i 

-0.0861* 

0.0315 

-2.7310 

R2 

F-  statistic 

0.1695 

4.5258* 

Adj-R2 

0.1321 

Diagnostic  tests 

Test 

F-statistic 

Prob.  value 

x2NORMAL 

0.1431 

0.8800 

x2serial 

0.2124 

0.8677 

X2ARCH 

0.0485 

0.8856 

X2  WHITE 

0.2996 

0.6014 

X2REMSAY 

2.9919 

0.1162 

Note:  x 2  NORMAL  is  for  normality  test,  x2  SERIAL  for  LM  serial  correlation  test,  x2 
ARCH  for  autoregressive  conditional  heteroskedasticity,  x2  WHITE  for  white 
heteroskedasticity  and  x2  REMSAY  for  Resay  Reset  test. 

*  The  significant  at  1%. 

**  The  significant  at  5%. 

»»»  j|ie  significant  at  10%  levels. 


Fig.  1.  Plot  of  cumulative  sum  of  recursive  residuals. 


Fig.  2.  Plot  of  cumulative  sum  of  squares  of  recursive  residuals. 


Table  7 

Chow  forecast  test 


Chow  forecast  test:  Forecast  from  1993Q1  to  2010Q4 

F-statistic  0.6843  Probability  0.9389 

Log  likelihood  ratio  82.889  Probability  0.1787 


capita  have  positive  and  negative  signs  (inverted-U  shaped  rela¬ 
tion)  on  C02  emissions  but  both  are  statistically  insignificant.  The 
impact  of  financial  development  is  positive  on  C02  emissions  and 
it  is  statistically  significant  at  10%  significance  level.  Electricity 
consumption  is  inversely  linked  with  C02  emissions  at  5%  level  of 
significance.  Trade  openness  increases  C02  emissions  at  5%  level  of 
significance.  The  statistically  significant  estimate  of  lagged  error 
term  i.e.  £CMt_!  with  negative  sign  corroborates  our  established 
long  run  relationship  between  carbon  emissions,  industrialization, 
electricity  consumption,  financial  development  and  trade  open¬ 
ness.  The  lagged  error-correction  term  is  negative  and  statistically 
significant.  This  implies  that  short  run  deviations  from  short  run 
towards  long  run  are  corrected  by  8.61%  in  each  quarter  and  it 
would  take  almost  four  years  to  reach  the  stable  log  run  equili¬ 
brium  path  in  C02  emissions  model  in  case  of  Bangladesh. 

The  results  of  diagnostic  tests  suggest  that  the  short  run  model 
passes  all  tests  successfully  such  as  autoregressive  conditional 
heteroskedasticity,  white  heteroskedasticity  and  specification  in 
the  model.  This  indicates  that  there  is  no  problem  of  autoregres¬ 
sive  conditional  heteroskedasticity.  The  variance  is  homoscedastic 
and  functional  form  of  the  short  run  model  is  well  organized.  The 
empirical  finding  shows  that  short  run  empirical  evidence  is 
consistent  and  stable  for  policy  purpose  regarding  carbon  emis¬ 
sions  in  case  of  Bangladesh. 

The  stability  of  the  ARDL  bounds  testing  estimates  is  investi¬ 
gated  by  applying  the  CUSUM  and  CUSUMsq  tests.  The  results  are 
shown  in  Figs.  1  and  2.  The  plots  of  the  CUSUM  statistics  are  well 
within  the  critical  bounds. 

The  plots  of  the  CUSUM  of  squares  statistics  are  not  well  within 
the  critical  bounds.  Furthermore,  we  apply  Chow  forecast  test  to 
examine  the  significance  of  structural  breaks  in  an  economy  for 
the  period  1993-2005  (explanations  are  given  above).  In  this 
study,  F-statistics  computed  in  Table  7  suggests  that  there  is  no 
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Table  8 

Variance  decomposition  method 


Period 

In  Ct 

In  It 

In  Et 

In  Ft 

In  Of 

Variance  decomposition  of  In  Ct 

1 

100.0000 

0.0000 

0.0000 

0.0000 

0.0000 

2 

99.5455 

0.1750 

0.0153 

0.2230 

0.0409 

3 

98.5077 

0.5697 

0.0227 

0.6972 

0.2025 

4 

96.7884 

1.1948 

0.0866 

1.3590 

0.5710 

5 

94.4321 

1.9645 

0.3087 

2.1327 

1.1617 

6 

91.5015 

2.7499 

0.8520 

2.9624 

1.9341 

7 

88.0636 

3.4399 

1.8918 

3.8088 

2.7956 

8 

84.1943 

3.9708 

3.5510 

4.6476 

3.6362 

9 

80.0078 

4.3296 

5.8356 

5.4633 

4.3635 

10 

75.6685 

4.5409 

8.6151 

6.2481 

4.9271 

11 

71.3684 

4.6474 

11.661 

7.0006 

5.3223 

12 

67.2851 

4.6936 

14.7218 

7.7233 

5.5760 

13 

63.5489 

4.7164 

17.5885 

8.4185 

5.7275 

14 

60.2324 

4.7414 

20.1275 

9.0845 

5.8140 

15 

57.3570 

4.7834 

22.2807 

9.7155 

5.8632 

16 

54.9070 

4.8482 

24.0492 

10.3025 

5.8929 

17 

52.8427 

4.9356 

25.4724 

10.8364 

5.9126 

18 

51.1109 

5.0412 

26.6103 

11.3108 

5.9266 

19 

49.6526 

5.1592 

27.5292 

11.7236 

5.9356 

20 

48.4100 

5.2837 

28.2924 

12.0744 

5.9392 

Variance  decomposition  of  In  It 


1 

0.0180 

99.9819 

0.0000 

0.0000 

0.0000 

2 

0.2323 

98.4507 

0.6372 

0.0764 

0.6031 

3 

0.5514 

94.4658 

2.3904 

0.9652 

1.6269 

4 

0.6882 

86.0900 

6.6907 

2.2462 

4.2846 

5 

0.6252 

74.1082 

13.4492 

3.3848 

8.4324 

6 

0.5005 

61.2797 

21.3213 

3.9109 

12.987 

7 

0.3998 

50.2847 

28.6677 

3.9119 

16.7357 

8 

0.3253 

42.0615 

34.7154 

3.6631 

19.2345 

9 

0.2784 

36.3654 

39.3841 

3.3667 

20.6052 

10 

0.2912 

32.5830 

42.8654 

3.1109 

21.1493 

11 

0.4078 

30.1260 

45.3965 

2.9195 

21.1499 

12 

0.6553 

28.5328 

47.1995 

2.7905 

20.8217 

13 

1.0288 

27.4722 

48.4722 

2.7133 

20.3133 

14 

1.4946 

26.7220 

49.3869 

2.6747 

19.7216 

15 

2.0047 

26.1430 

50.0860 

2.6600 

19.1060 

16 

2.5127 

25.6543 

50.6783 

2.6551 

18.4994 

17 

2.9840 

25.2126 

51.2383 

2.6478 

17.9171 

18 

3.4002 

24.7968 

51.8099 

2.6292 

17.3637 

19 

3.7565 

24.3984 

52.4120 

2.5949 

16.8380 

20 

4.0581 

24.0154 

53.0461 

2.5439 

16.3363 

Variance  decomposition  of  In  Et 

1 

3.9586 

1.9912 

94.0501 

0.0000 

0.0000 

2 

3.2369 

2.4101 

94.3385 

0.0093 

0.0049 

3 

2.9013 

2.2918 

94.7052 

0.0959 

0.0055 

4 

2.3851 

2.3525 

94.9682 

0.2639 

0.0300 

5 

1.8624 

2.5721 

94.9398 

0.5394 

0.0860 

6 

1.5526 

2.9444 

94.4067 

0.9161 

0.1799 

7 

1.6580 

3.4358 

93.2158 

1.3789 

0.3112 

8 

2.2810 

4.0012 

91.3460 

1.8994 

0.4721 

9 

3.3867 

4.5902 

88.9289 

2.4456 

0.6484 

10 

4.8272 

5.1588 

86.2003 

2.9897 

0.8238 

11 

6.4054 

5.6780 

83.4190 

3.5127 

0.9848 

12 

7.9399 

6.1340 

80.7985 

4.0046 

1.1227 

13 

9.3031 

6.5251 

78.4754 

4.4617 

1.2345 

14 

10.4293 

6.8558 

76.5106 

4.8833 

1.3209 

15 

11.3044 

7.1336 

74.9076 

5.2696 

1.3846 

16 

11.9487 

7.3667 

73.6340 

5.6209 

1.4295 

17 

12.4003 

7.5631 

72.6390 

5.9377 

1.4597 

18 

12.7035 

7.7303 

71.8665 

6.2207 

1.4788 

19 

12.9004 

7.8753 

71.2629 

6.4714 

1.4898 

20 

13.0274 

8.0039 

70.7816 

6.6917 

1.4952 

Variance  decomposition  of  In  Ft 

1 

0.0866 

20.8889 

0.2791 

78.7453 

0.0000 

2 

0.7384 

15.8874 

0.0949 

83.2738 

0.0053 

3 

2.0547 

12.4193 

0.1382 

85.3799 

0.0077 

4 

3.5066 

9.9519 

0.4568 

86.0630 

0.0214 

5 

4.8811 

8.3185 

1.0580 

85.7045 

0.0376 

6 

6.0735 

7.2921 

1.9280 

84.6609 

0.0453 

M.  Shahbaz  et  al.  /  Renewable  and  Sustainable  Energy  Reviews  31  (2014)  575-586 


583 


Table  8  ( continued ) 


Period 

In  C, 

In  It 

In  Et 

In  Ft 

In  Of 

7 

7.0740 

6.7030 

3.0202 

83.1610 

0.0414 

8 

7.9027 

6.4277 

4.2844 

81.3479 

0.0370 

9 

8.5835 

6.3730 

5.6753 

79.3184 

0.0495 

10 

9.1346 

6.4627 

7.1581 

77.1521 

0.0922 

11 

9.5692 

6.6334 

8.7075 

74.9220 

0.1677 

12 

9.8987 

6.8362 

10.3038 

72.6926 

0.2685 

13 

10.1361 

7.0379 

11.9297 

70.5148 

0.3812 

14 

10.2967 

7.2200 

13.5670 

68.4239 

0.4922 

15 

10.3976 

7.3753 

15.1961 

66.4397 

0.5911 

16 

10.4563 

7.5043 

16.7956 

64.5708 

0.6728 

17 

10.4891 

7.6117 

18.3449 

62.8181 

0.7360 

18 

10.5102 

7.7035 

19.8253 

61.1782 

0.7826 

19 

10.5314 

7.7854 

21.2221 

59.6452 

0.8155 

20 

10.5613 

7.8623 

22.5260 

58.2122 

0.8379 

Variance  decomposition  of  In  Ot 

1 

2.3399 

9.5410 

0.0515 

1.0278 

87.0397 

2 

2.0852 

5.4774 

0.0763 

1.6057 

90.7551 

3 

2.0315 

3.1461 

0.3744 

1.8523 

92.5955 

4 

2.0556 

2.5843 

0.7437 

2.4767 

92.1394 

5 

2.0776 

3.2188 

1.1776 

3.5712 

89.9545 

6 

2.0449 

4.3972 

1.6508 

5.2813 

86.6256 

7 

1.9355 

5.5794 

2.1820 

7.5958 

82.7070 

8 

1.7676 

6.4726 

2.7940 

10.3633 

78.6023 

9 

1.5856 

6.9991 

3.5056 

13.3207 

74.5888 

10 

1.4340 

7.2137 

4.3181 

16.1877 

70.8463 

11 

1.3372 

7.2182 

5.2126 

18.7454 

67.4863 

12 

1.2935 

7.1098 

6.1538 

20.8756 

64.5670 

13 

1.2835 

6.9606 

7.0994 

22.5548 

62.1014 

14 

1.2836 

6.8158 

8.0103 

23.8232 

60.0668 

15 

1.2761 

6.6999 

8.8572 

24.7504 

58.4161 

16 

1.2536 

6.6231 

9.6237 

25.4103 

57.0891 

17 

1.2181 

6.5870 

10.3050 

25.8673 

56.0224 

18 

1.1770 

6.5880 

10.9061 

26.1721 

55.1566 

19 

1.1392 

6.6199 

11.4380 

26.3622 

54.4405 

20 

1.1115 

6.6760 

11.9150 

26.4647 

53.8326 

significant  structural  break  in  the  economy  during  the  sample 
period.  The  chow  forecast  test  is  more  reliable  and  preferable  than 
graphs  [88],  This  confirms  that  the  ARDL  estimates  are  reliable  and 
efficient. 

The  direction  of  causality  between  industrialization,  electricity 
consumption,  financial  development,  trade  openness  and  C02  emis¬ 
sions  was  found  by  applying  innovative  accounting  approach  (IAA) 
rather  than  the  VECM  Granger  causality  method.  The  VECM  Granger 
causality  is  suitable  to  detect  a  causal  relationship  between  the 
variables  within  the  sampled  period.  To  determine  causality  ahead 
the  sample  period,  the  innovative  accounting  approach  is  much  better. 
The  innovative  accounting  approach  is  combination  of  variance 
decomposition  and  impulse  response  function.  The  variance  decom¬ 
position  approach  indicates  the  magnitude  of  the  predicted  error 
variance  for  a  series  accounted  for  by  innovations  from  each  of  the 
independent  variable  over  different  time-horizons  beyond  the  selected 
time  period.  It  is  pointed  by  Pesaran  and  Shin  [79]  that  generalized 
forecast  error  variance  decomposition  method  shows  the  proportional 
contribution  in  one  variable  due  to  innovative  shocks  stemming  in 
other  variables.  The  main  advantage  of  this  approach  is  that  like 
orthogonalized  forecast  error  variance  decomposition  approach;  it  is 
insensitive  with  ordering  of  the  variables  because  ordering  of  the 
variables  is  uniquely  determined  by  VAR  system.  Further,  the  general¬ 
ized  forecast  error  variance  decomposition  approach  estimates  the 
simultaneous  shock  effects.  Engle  and  Granger  [89]  and  Ibrahim  [90] 
argued  that  with  VAR  framework,  variance  decomposition  approach 
produces  better  results  as  compared  to  other  traditional  approaches. 
The  results  of  variance  decomposition  approach  are  described  in 
Table  8.  The  empirical  evidence  indicates  that  a  48.41%  portion  of 
C02  emissions  is  contributed  by  its  own  innovative  shocks  and  one 


standard  deviation  shock  in  industrialization  explains  energy  pollu¬ 
tants  by  5.28%  which  is  minimal. 

Electricity  consumption  contributes  to  C02  emissions  by  28.29% 
due  to  one  standard  shock  stemming  in  energy  utilization.  The 
share  of  financial  development  and  trade  openness  in  C02  emis¬ 
sions  is  minimal  i.e.  12.07%  and  5.93%  respectively.  Electricity 
consumption  explains  industrial  development  53.04%  due  to 
innovative  shocks  in  electricity  consumption.  A  24.01%  of  indus¬ 
trialization  is  explained  by  own  standard  shock.  The  contribution 
of  C02  emissions,  financial  development  and  trade  openness  to 
industrialization  is  4.05%,  2.54%  and  16.33%  due  to  one  standard 
shock  arises  in  these  series  respectively.  Electricity  consumption  is 
contributed  13.02%,  8.00%,  6.16%  and  1.49%  by  C02  emissions, 
industrialization,  financial  development  and  trade  openness  while 
rest  i.e.  70.78%  is  contributed  by  its  own  sakes. 

Similarly,  a  major  contribution  (58.21%)  in  financial  development 
by  its  own  standard  shock  while  contribution  by  C02  emissions, 
industrialization  and  trade  openness  is  10.56%,  7.86%  and  0.83% 
respectively.  The  contribution  in  financial  development  is  due  to 
one  standard  shock  in  electricity  consumption  is  22.52%.  C02  emis¬ 
sions,  industrialization,  electricity  consumption  contribute  to  trade 
openness  by  1.13%,  6.61%  and  11.53%  respectively.  Financial  develop¬ 
ment  explains  trade  openness  by  26.36%  through  one  standard 
deviation  shocks  stems  in  financial  development.  Overall  our  results 
indicate  that  electricity  consumption  causes  C02  emissions.  Indus¬ 
trialization  is  also  a  cause  of  electricity  consumption.  The  unidirec¬ 
tional  causality  is  running  from  electricity  consumption  to  financial 
development  and  financial  development  causes  trade  openness. 

The  impulse  response  function  (Fig.  3)  is  alternative  to  variance 
decomposition  method  shows  how  long  and  to  what  extent 
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Fig.  3.  Impulse  response  function. 
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dependent  variable  reacts  to  shock  stemming  in  the  independent 
variables.  The  results  indicate  that  the  response  in  C02  emissions  due 
to  forecast  error  stemming  in  industrialization  initially  rises,  goes  to 
peak  and  then  starts  to  decline  after  8th  time  horizon.  This  presents 
the  phenomenon  of  environmental  Kuznets  curve  between  industria¬ 
lization  and  C02  emissions.  The  contribution  of  financial  development 
to  C02  emissions  is  positive  till  10th  time  horizon.  C02  emissions 
respond  negatively  but  positively  after  5th  time  horizon  due  to 


forecast  error  in  electricity  consumption.  The  response  in  C02  emis¬ 
sions  is  negative  after  4th  time  horizon  due  to  trade  openness. 
Electricity  consumption  and  trade  openness  contribute  positively  to 
industrialization. 

The  response  in  electricity  consumption  is  positive  due  to 
forecast  error  arising  in  industrial  and  financial  development. 
Financial  development  responds  positively  due  to  shocks  in  C02 
emissions,  industrialization  and  electricity  consumption  but 
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negatively  due  to  shocks  occurring  in  trade  openness.  The 
response  in  trade  openness  is  inverted  U-shaped  due  to  shocks 
in  C02  emissions  and  same  for  trade  openness  and  industrializa¬ 
tion  but  positive  due  to  forecast  error  in  electricity  consumption. 
Financial  development  inversely  contributes  to  trade  openness. 

5.  Conclusion  and  policy  implications 

This  paper  deals  with  empirical  investigation  of  the  relation¬ 
ship  between  C02  emissions  and  industrialization  by  incorporat¬ 
ing  electricity  consumption,  financial  development  and  trade 
openness  as  potential  determinants  of  C02  emissions  in  case  of 
Bangladesh.  We  have  applied  structural  break  unit  root  test  and 
the  long  run  relationship  between  the  variables  is  investigated  by 
applying  the  ARDL  bounds  testing  approach  to  cointegration  in  the 
presence  of  structural  break  in  the  series.  Causal  relationship 
between  C02  emissions,  industrialization,  electricity  consumption, 
financial  development  and  trade  openness  is  scrutinized  by  using 
innovative  accounting  approach. 

Our  results  find  that  the  environmental  Kuznets  curve  exits 
between  industrialization  and  C02  emissions  in  Bangladesh.  Electri¬ 
city  consumption  adds  in  C02  emissions.  Financial  development 
increases  carbon  emissions.  Trade  openness  and  energy  pollutants 
are  positively  linked.  The  causal  analysis  reveals  that  electricity 
consumption  causes  C02  emissions  and  industrialization.  The  find¬ 
ings  of  this  paper  have  an  important  attention  in  Bangladesh  where 
the  country  facing  a  serious  climate  change  problem.  Based  on  the 
findings,  this  paper  suggests  that  Bangladesh  government  should 
need  to  take  necessary  policy  to  improve  the  environmental  quality. 
In  order  to  improve  the  environmental  quality,  government  should 
consider  the  sustainable  energy  policy  that  will  accelerate  the 
renewable  energy  production  and  reduce  the  energy  subsidy  in  the 
price  level  that  will  mitigate  C02  emissions.  Government  should 
more  focus  on  the  clean-energy  investment/green  investment  by 
adding  the  implemention  of  environmental  regulation  for  the  dirty 
industries.  In  addition  regulatory  institution  provides  the  financial 
incentive  for  clean  technology.  Financial  development  causes  trade 
openness  and  electricity  consumption.  Hence,  sound  and  developed 
financial  system  that  can  attract  investors,  boost  the  stock  market 
and  improve  the  efficiency  of  economic  activities  should  be  encour¬ 
aged  in  the  country.  Sound  financial  sector  also  improves  the  overall 
trade  situation  of  the  country  because  the  cost  of  doing  business  is 
low.  Dirty  industrialization  leads  to  the  increase  of  the  emissions  in 
the  country.  In  Bangladesh,  leather,  chemical  and  shipbuilding 
industries  use  large  amounts  of  toxic  chemicals  that  has  serious 
negative  consequences  on  the  environment.  Government  should 
regulate  these  industries  and  advise  them  to  adopt  environment 
friendly  technology  to  enhance  production  of  leather,  chemical  and 
shipbuilding. 
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